DB 755 PCT 



60 



CLAIMS 

1. Equalization method for the cancellation of isofrequential interferers in 
signals received by array antennas (ARY) of base stations of a cellular mobile system, 
said signals (y(w)) being received under the form of bursts of data symbols of a 

5 transmission burst, including a training sequence ({*„}) known to the receiver and 

usable to estimate the impulse response of the transmission channel, said bursts 
received being demodulated, digitized, stored, and synchronized through discrimination 
of the training sequence, including the following steps: . 

a) spatial filtering of signals coming from the antennas (yM) » using a spatial filter 
10 (BMF) that combines among them corresponding n-th symbols of the bursts on the 

basis of n-th weights of a sequence of weights (w opf (n)) obtaining a beamformed 
received sequence (y^n)); 

b) time filtering (CAMP) data of the training sequence ({*„}) for estimating a 
sequence of samples (h(n)) of the impulse response of the transmission channel; 

15 c) calculating a square error (e n ) originated from the comparison between the 
corresponding n-th samples of the estimate impulse response (h(w)) and the 
beamformed sequence (y^C*)); 

d) joint optimization of said weights and samples (w^h), h opt (n)) by minimizing a 
functional (FUNZ) corresponding to said mean square error (e„); 
20 e) time filtering of the signal (y^W) generated by said spatial filter (BMF) using a 
filter (FADAT) matched to the sequence of samples (h opf (w)) of the impulse 

response of the transmission channel; 

f) providing said sequence of samples of the impulse response of the transmission 

channel {h opt {nij and a sequence of sample at the output of said matched filter 
25 (FADAT) to a Viterbi estimator for the detection of the received data; 

g) updating the optimized values of said weights and samples (yv opt {n), h o;7 ,( n )). out 
of the time slot of the training sequence, through repetition of the previous steps 
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a)...f), replacing the symbols of the training sequence ({*„}) by a low delay 

estimated sequence ({-?'„ }) , thus rendering said equalization adaptive, 

characterized in that while the previous steps are executed a parallel continuous 
monitoring of the power received on the antennas (ARY) for detecting possible peaks 

5 exceeding a given threshold in coincidence of the arrival of a strong non-stationary 
interferer is carried out, in case of detection of a peak over threshold at an instant n t 
step d) is modified for the all duration starting from the instant n up to an instant 
n + AT of a transient devoted to the rejection of strong interferent outside the 
midamble, by updating the sole weights (w op ,(*)) of the spatial filter (BMF) maintaining 

10 the samples of the impulse response (h^(«)) constant at the values they have at the 
instant n , then the execution of unmodified step d) is restored at the end of said 
transient. 

2. Equalization method according to claim 1 , characterized in that in absence 
of a strong non-stationary interferer the preceding step d) is executed in closed form by 
15 the following sub-steps: 

calculation of a spatial covariance matrix R^, of the signals coming from the 

antennas (y(w)), a covariance matrix of the data samples ({x n }, {x' n }) , a 
cross-covariance matrix R^ between data samples and inputs to the antennas 
({jc n },y(w)), and the Schur complement matrix R, 1 = R^ - R^R^R^ , being * the 



20 complex conjugate, of a matrix R = 



R yy R W 



estimation of an eigenvector h opt (n) of matrix R^ associated to the minimum 
module eigenvalue y given by h*Rj*h = -?\ said eigenvector including said 
optimized samples; 

estimation of said optimized weights by calculating a vector w opr (*) solution of 

25 the equation: R^w^n) =R jy h tfp , („). 

3. Equalization method according to claim 1, characterized in that in 
presence of a strong non-stationary interferer the preceding step d) is executed in 
closed form by the following sub-steps: 




calculation of a spatial covariance matrix R^(« + AT) of the signals coming 
from the antennas (y{n)) and a complex conjugate cross-covariance matrix 
R^yin) between data samples and inputs to the antennas ({* n },y(«)), matrix 
+ AT) being calculated at instant n , and then kept constant for the whole 
AT , using vectors y{n) obtained from the samples of signals on the antennas 
present in the stored bursts, starting from instant n up to the instant n-t- AT, while 
R^(rt) being calculated at each instant n-th of the transient; 

estimation of said optimized weights by calculating a vector w op ,(w) solution of 
the equation: R„(n + AT) w = *C - where h opt is a vector including 

the samples of the impulse response at the instant n . 

4. Equalization method according to claim 1 , characterized in that in absence 
of a strong non-stationary interferer the preceding step d) is executed in iterative form 
by the following sub-steps: 
d1) assignment of arbitrary values, preferably null, to the elements of a structured 

'L(w) 0 



triangular matrix L(w), defines as follows: L(w) = 



U(/i) L»_ 



, wherein L(/z) 



and L,(rt) are triangular matrices and U(n) is a square matrix, or more in 

general a rectangular one; 
d2) assignment of arbitrary values, one of which non null, to the elements of a vector 
h^(«) storing said optimized samples of the impulse response of the 

transmission channel; 
d3) assignment of arbitrary values, one of which at least non null, to the elements of 
a vector w opf (n) storing said weights of the spatial filter (BMF); 

d4) updating of said structured matrix E(n) applying at instant n-th the following 

iterative relation: \L(n) o] = [E(n-l) g„]Q("). where g„ is the following 



updating vector g n = 



structured in two component vectors, out of which: 



y(") = [^i(").^2(") yM( n )] T ' ls a f,rst vector whose M elements are symbols 

coming from M antennas at instant n -th, present in a same number of stored 




bursts, and x(n)^[x n (n) J x„^(n\.... J x n ^ L ^(n)] T is a second vector whose L 
elements are stored symbols of the transmitted training or estimated sequence 
{*'„}) respectively, assuming a duration of the impulse response equal to 
L symbols; being Q(n) a matrix transformation of unitary rotation, implemented 
through a sequence of elementary transformations Q fl ,Q I} that applied to 

[l(az-1) g n ] produces the cancellation of column g n introducing a zero for 

each elementary transformation; 
d5) determination of a weight vector h n through the following expression: 

L s (n)V s (n)h n = h op/ (*-l). according to which said vector of optimum 



10 samples is updated as follows h («) = 



d6) updating of said optimum weight vector through the following expression: 
L'Ww^W = U»h c ». 

5. Equalization method according to claim 4, characterized in that said 
expression calculated at step h5) is divided into two calculation phases, introducing an 

15 unknown intermediate vector u n defined by the following first equation: 

L*(>*)ti n = u„, being u rt a solution of the following second equation: 
L x (n)u n = h opt (n-l) that replaced in its turn in the starting equation gives h n , 

simplifying the relevant calculation. 

6. Equalization method according to claim 4, characterized in that said matrix 

20 L(n - 1) is multiplied by VX , where X is an arbitrary real constant, lower than one 
unit, selected in relation to the non stationariness of the channel. 

7. Adaptive equalization method according to claim 1 , characterized in that, 
in presence of a strong non-stationary interferer the preceding step d) is executed in 
iterative form by the following sub-steps: 

25 hi) advanced updating within the symbol time n of a sole triangular matrix L 
belonging to a structured triangular matrix L(n), defined as follows: 

"L(n + Ar) 0 
U(n) L,(i.)J" 

and U{n) is a square matrix, or more in general a rectangular one; the advanced 



L(„)= 



.wherein L»(« + Ar) and L^n) are triangular matrices 




updating of said matrix L being executed using vectors y(n) obtained from the 
samples of signals on the antennas, present in the stored bursts, starting from 
instant n up to an instant n + AT\ matrix L being calculated at instant n and 
then kept constant for the whole AT ; 
h2) starting from instant n up to an instant n + AT , updating of matrices U(w) and 
L,(«) of said triangular matrix L(n) applying at instant n-th the following 
iterative relation: \L(n) o] = [E(n-l) g„]Q(w), where g„ is the following 



updating vector g n = 



structured in two component vectors, out of which: 



y(n) = [yXn\y 2 (n) 9 ....,y M (n)] T \s a first vector whose M elements are symbols 
coming from M antennas at instant n -th, present in a same number of stored 
bursts, and x(n) = [^ /I (/7),^ n+1 (n),....,x n+I . l (n)] r is a second vector whose L 
elements belong to said low delay estimated sequence ({*'„}), assuming a 

duration of the impulse response equal to L symbols; being Q(n) a matrix 
transformation of unitary rotation, implemented through a sequence of elementary 
transformations Q 0 ,Q|.--, that applied to [L(w-l) g n ] produces the 
cancellation of column g„ introducing a zero for each elementary transformation; 
h3) updating of said optimum weight vector through the following expression: 
U(n +AT)v( opt {n) = \J*(n)h opt , where h opt is a vector including the samples of 

the impulse response at the instant n . 
8. Equalization method according to claim 7, characterized in that said 
advanced updating of the sole triangular matrix L , is made applying the generic said 
elementary transformation Q, to the first M elements of the /-th column of said 

triangular matrix L(w), where M is the order of matrix L, repeating the elementary 
transformation Q, for all the values of index / and storing the transformations occurred. 
<J 9. Equalization method according to claim 7^p-8, characterized in that said 
matrix L(« - l) is multiplied by VX , where A is an arbitrary real constant, lower than 
one unit, selected in relation to the non stationariness of the channel. 

10. Equalization method according to — — tho prov i ouo ■■ c t a i rno , 
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characterized in that, coinciding with any new symbol of low delay estimated sequence 
({£•„})., the L elements of said vector x(«) are replaced by a same number of symbols 
copied as a whole from the path selected in the trellis at the present symbol time, 
preventing by this the propagation of errors in the joint estimate of said optimum 

5 weights and channel impulse response samples (w 

opt 

a 11. Equalization method according to^any^ nf- ahe — pfeee^awp-wtotms, 

characterized in that when said training sequence (f*?j) is included in a midamble of 

the transmission burst, said equalization of the step g) starts from a first end of the 
midamble and proceeds towards the end of the burst, previously stored and 

10 synchronized, and afterwards it restarts from the second end of the midamble towards 
the other end of the burst, using estimated data of the first semiburst to have a longer 
group of estimated data for the equalization of the second semiburst. 

12. Equalizer in receivers of base stations of a cellular mobile system, in 
particular for the cancellation of isofrequential interferers in signals received by array 

15 antennas (ARY) of base stations of a cellular mobile system, said signals (y(n)) being 
received under the form of bursts of data symbols of a transmission burst, including a 
training sequence ({*„}) known to the receiver and usable to estimate the impulse 
response of the transmission channel, said bursts received being demodulated, 
digitized, stored, and synchronized through discrimination of the training sequence, 

20 including: 

- a spatial filter (BMF) of signals coming from the antennas (y(w)) that combines 
among them corresponding n-th symbols of the bursts on the basis of n-th weights 
of a sequence of weights (w op ,(n)) obtaining a beamformed received sequence 

(y «*(")): 

25 - a time filter (CAMP) for filtering data of the training sequence ({x n }) and estimating 
a sequence of samples (h(w))of the impulse response of the transmission 
channel; 

- means (z) suitable for calculating a square error (e„) originating from the 
comparison between the corresponding n-th samples of the estimate impulse 




response (h(«))and the beamformed sequence (y ouf (n)) ; 

- means (FUNZ) suitable for the joint optimization of said weights and samples 
(v/ opl {n), h opt {nj) and to minimize a functional (FUNZ) corresponding to said mean 
square error (e„); 

5 - a time matched filter (FADAT) of the signal (y oa («)) generated by said spatial filter 
(BMF). matched to the sequence of samples (h opt (nj) of the impulse response of 
the transmission channel; 

- a Viterbi estimator (STIM) for estimating the symbols ({*„}) of the transmission 
burst corresponding to the sequence coming out from the time matched filter 

10 (FADAT); 

- means (DECTEMP) suitable for replacing out of the time slot of the training 
sequence the symbols of the training sequence ({*„}) with a low delay estimated 

sequence ({*'„}) . thus rendering said equalization adaptive, 

characterized in that it further includes: 
15 - means (PICRIV) suitable for the continuous monitoring of the power received on 
the antennas (ARY) and for detecting possible peaks exceeding a given threshold 
in coincidence of the arrival of a strong non-stationary interferer; 

- selection means (PICSEL, SWT) controlled by said continuous monitoring means 
(PICRIV) suitable for the selection of said means (SYSTOL1, 3, 4, 5) for the joint 

20 optimization of weights and samples, in case a detection of a peak over threshold 

doesn't occur, and for the selection of means (SYSTOL2, SYSTOL5) devoted to 
the rejection of a strong interfered in case of detection of a peak over threshold 
occur at an instant n outside the midamble; the latter means being operative for 
the all duration starting from the instant n up to an instant n + AT of a rejecting 
transient and being arranged for updating the sole weights (w op ,(n)) of the spatial 
filter (BMF) maintaining the samples of the impulse response (h op ,(w)) constant at 
the values they have at the instant n; the preceding means for the joint 
optimization (SYSTOL1, 3, 4, 5) being restored at the end of said rejection 
transient. 

30 13. Equalizer according to claim 12, characterized in that in absence of strong 



25 
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10 L(/j) = 



image of said first systolic array (SYSTOL1) in which the 



non-stationary interferer said joint optimization means (SYSTOL1, 3, 4, 5) are 
constitute by the following means: 

- first delay means (RIT) for delaying with a fixed delay (D) an algebraic vector y{n) 
of M elements, as many as the antennas, whose elements are symbols of the 
stored bursts of the signal coming from the antennas; 

- a first systolic array (SYSTOL1) receiving said vector y(n) and symbols of said 

training sequence ({*„}) or estimated low delay ones, organized as in an 

algebraic vector x(n) of L elements, equal to the length of the channel impulse 
response, and updating, coinciding with each symbol time n-th, a triangular matrix 
'L(n) 0 

elements L(/i) and L,(/i) are also triangular matrices, and U(n) is rectangular, 
or square; 

- a second systolic array (AUTOVET), said second systolic array receiving said 
element L x (n) and updating a vector h^in) whose elements correspond to the 

15 optimized samples of the impulse response, 

- a third systolic array (SYSTOL5), placed in cascade to the second, receiving said 
elements h{n) and U(n), the vector h^in), and updating a vector w opt {n) 

whose elements correspond to the optimized weights. 

14. Equalizer according to claim 13, characterized in that said second 
20 systolic array (AUTOVET), is further divided in the cascade of two systolic arrays 
(SYSTOL3, SYSTOL4) having identical triangular structure, of which the upstream 
array (SYSTOL3) receives said element L f W, and the previous vector h opt {n-l) of 
the downstream array, and calculates an intermediate vector u n supplying it to the 
downstream array (SYSTOL4), reached also by element L*(w), for the calculation of 
25 said vector h opt (n) . 

15. Equalizer according to claim 12, characterized in that in presence of strong 
non-stationary interferer said joint optimization means (SYSTOL2, SYSTOL5) are 
constitute by the following means: 

- first delay means (RIT) for delaying with a fixed delay (D) an algebraic vector y{n) 
30 of M elements, as many as the antennas, whose elements are symbols' of the 
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stored bursts of the signal coming from the antennas; 
- a first systolic array (SYSTOL2) receiving said vector y(«) and the symbols of said 

low delay estimated sequence ({*'„}) organized as in an algebraic vector x{n) of 
L elements, equal to the length of the channel impulse response, and updating, 



image of said first systolic array (SYSTOL2) in which the elements L(n) and L> s {n) 
are also triangular matrices, and U(w) is rectangular or square matrix ; 
- second delay means (FIFO) suitable for delay element L calculated at instant n 
for the whole AT ; 

10 - a second systolic array (SYSTOL5), placed in cascade to the first, receiving said 
elements L(n) and tl(/z) , a vector h opt including the samples of the impulse 
response at the instant n , and updating a vector w opf (w) whose elements 

correspond to the optimized weights 
O 16. Equalizer according to claim 13 -o r 15 , characterized in that said vectors 

15 y(«) and x(n) reach in column, in the order indicated, the input of the cells of said first 
systolic array (SYSTOL1 , SYSTOL2) corresponding to the elements of the first column 
of said matrix h(n) , the whole of cells of the triangular matrix including the calculation 
means of the following matrix relation: 



20 where: 

- a pointed by the arrow indicates the new content of the cell; 

- A is an arbitrary real constant <, 1, selected in relation to the non stationariness of 
the channel. 

- oaf is an output signal; 

25 - in is a signal entering each cell; and 

- p is a both input and output signal, calculated as follows: 



5 



coinciding with each symbol time n-th, a matrix h(n) = 



'h{n + AT) 0 
U(«) L» 




in 




DB 755 PCT 



69 



# 



containing the alternative calculation means, such for which the signal p is only an 
output one, the signal out is absent, and the content a is updated as follows: 



